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We study gluon scattering amplitudes/ Wilson loops in Af = 4 super Yang-Mills 
| theory at strong coupling by calculating the area of the minimal surfaces in AdS^ 

based on the associated thermodynamic Bethe ansatz system. The remainder func- 
^ ■ tion of the amplitudes is computed by evaluating the free energy, the T- and Y- 

' functions of the homogeneous sine-Gordon model. Using conformal field theory 

(CFT) perturbation, we examine the mass corrections to the free energy around 
the CFT point corresponding to the regular polygonal Wilson loop. Based on the 
relation between the T-functions and the (/-functions, which measure the boundary 
entropy, we calculate corrections to the T- and Y-functions as well as express them 
at the CFT point by the modular S-matrix. We evaluate the remainder function 
around the CFT point for 8 and 10-point amplitudes explicitly and compare these 
analytic expressions with the 2-loop formulas. The two rescaled remainder functions 
show very similar power series structures. 
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1. Introduction 

It has been recognized that there exists an integrability structure in gluon scattering 
amplitudes in M = 4 super Yang-Mills theory at strong coupling. Gluon scattering 
amplitudes are dual to the Wilson loops made of light-like segments [HE]. By us- 
ing the AdS-CFT correspondence, the amplitudes at strong coupling are shown to 
be equal to the area of the minimal surface in AdS with the same null polygonal 
boundary [TJ. For n(> 6)-point amplitudes [3j[I] they are shown to differ from the 
Bern-Dixon-Smirnov (BDS) conjecture |5j. The deviation from the BDS conjecture 
is called the remainder function, which is a function of the cross-ratios of the external 
gluon momenta. 

It has been found that the equations for the minimal surface reduce to the Hitchin 
system and the area of the surface is determined by the Stokes data of its solutions 
[OHH] • The Stokes data and their cross-ratios obey the functional equations called the 
T-system [TU] and the Y-system [TT]. The area of the minimal surface is determined 
by the thermodynamic Bethe ansatz (TBA) equations [12] associated with the Y- 
system [8]. 



Y-systems are closely related to two-dimensional integrable models. For the 6- 
point gluon scattering amplitude, which corresponds to the minimal surface with 
a hexagonal null boundary in AdS^, the related Y-system and the TBA equations 
are those of the Z4-symmetric integrable model [7]. In our previous paper [13], we 
solved the TBA equations with chemical potential by the integrable perturbation 
of conformal field theory (CFT) and evaluated the free energy. In order to obtain 
the analytic form of the remainder function near the CFT point corresponding to 
the regular polygonal Wilson loop, it is moreover necessary to calculate small mass 
expansion of the Y-functions, which was determined numerically in [IB]. 

A key observation for computing the Y-functions analytically is as follows: The 
TBA free energy is obtained from the partition function on a cylinder with the 
periodic boundary condition. We can also consider the free energy on a cylinder with 
different boundary conditions. Affleck and Ludwig [2] introduced the (^-functions 
from this free energy, which count the ground state degeneracy of the system with 
boundaries (the boundary entropy). In [EKE], Dorey et al. studied the exact off- 
critical (^-functions for the purely elastic scattering theory and derived the integral 
equation for them. They evaluated the g-functions using the integrable perturbation 
of the boundary CFT. Remarkably, ratios of the (^-functions obey the same integral 
equations for the T-functions. This enables us to obtain the analytic expression 



1 



for the Y-functions from the (/-functions and determine the analytic form of the 
remainder function near the CFT point. 

In this paper, we study the remainder function for 2n-point gluon scattering 
amplitudes at strong coupling, which correspond to minimal surfaces in Ad S3 with a 
2n-gonal light-like boundary. This corresponds to the case where the gluon momenta 
are in R 1 ' 1 . In [9] the related integrable system was shown to be the homogeneous 
sine-Gordon model (HSG) [17] with purely imaginary resonance parameters. The 
relevant CFT is the generalized parafermions SU(n — 2)2/U(l) n ~ 2 [18]. In this paper 
we will study the boundary and bulk perturbation of the generalized parafermions 
and calculate the ratios of the (/-functions. For the octagon (n = 4) and the decagon 
(n = 5), we will calculate the perturbative corrections to the T-/ Y-functions, the free 
energy and the remainder function explicitly. We compare these analytic expressions 
of the remainder function with the 2-loop formulas proposed in [T5H2"Tj around the 
CFT limit. 

The above analytic results are important to understand the structure and the 
momentum dependence of the amplitudes at strong coupling exactly. The purpose 
of this paper is to take a step toward this direction by analyzing the TBA system 
near the CFT point from the conformal perturbation theory (CPT). A point in our 
discussion is that not only the free energy but also the Y-functions can be discussed in 
this framework owing to the relation between the T-/Y-functions and the (/-functions. 

This paper is organized as follows. In section 2, we review the TBA system for 
the minimal surfaces in AdS%. In section 3, we discuss the HSG model and the 
free energy for the integrable bulk perturbations. In section 4, we study the HSG 
model from the CPT of the generalized parafermions and the relation between the T- 
functions and the g- functions. In section 5, we investigate the small mass expansions 
of the (/-function and the remainder function for the octagon. In section 6, we argue 
the corrections to the free energy, the T-/ Y-functions and the remainder function 
for the decagon. In section 7, we compare the analytic expressions for the remainder 
functions for the octagon and the decagon with the 2-loop results. In section 8, we 
give a summary and a discussion. In Appendix A, we present the high-temperature 
expansion of the (/-function for the Ising model in detail. In Appendix B, we study 
the CPT of an SU(2) coset model with fractional level. In Appendix C, we discuss 
the expansion of the Y-functions in the case of complex masses. In Appendix D, we 
examine the structure of the expansion of the remainder functions at higher orders. 
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2. Review of TBA system for minimal surfaces in AdS 3 

2.1. Problem 

Gluon scattering amplitudes at strong coupling are evaluated by the area of minimal 
surfaces in AdS§ [llj. The minimal surfaces have a polygonal boundary located on the 
AdS boundary. Each edge of the polygonal boundary represents a gluon momentum. 
Throughout this paper we consider the case where the minimal surfaces stretch inside 
the maximal AdS? subspace. In this case, the polygonal boundary is contained 
inside R 1 ' 1 at the AdS boundary. Namely, we consider amplitudes with all the gluon 
momenta being restricted in M 1 ' 1 . For general 2n-point amplitudes, we label the 
vertices of the polygon in the light-cone coordinates as x 2 k~i = ( x t^ x k)^ x zk = 
( x t+ii x k)i k = l)---)^ (with identification x h+1 = xf). The gluon momenta are 
then given by Xj — Xj+i, j = 1, . . . 2n (x2n+i — x\). These data completely fix the 
shape and the area of the minimal surfaces. 

As the gluon scattering amplitudes in M = 4 SYM are infrared divergent, the area 
of the minimal surfaces are also divergent. Nevertheless, the structure of the diver- 
gence has been well studied and one can read off meaningful results after specifying 
an appropriate regularization scheme. The finite results are commonly analyzed in 
the form of the remainder function, namely the deviation of the area from the con- 
jecture of Bern, Dixon and Smirnov j5]. It is a dual-conformally invariant quantity 
and hence a function of cross-ratios of Xj. The form of the function is of our central 
interest. 

2.2. T-/Y-system and TBA equations 

Despite the difficulty in analytically constructing the minimal surfaces for general 
polygonal boundaries, it turned out that the area can be computed by solving TBA 
type integral equations. This subsection is devoted to a brief summary of the TBA 
equations and the associated T-/Y-system for the minimal surfaces in AdS?,. For 
details, see 0[El[9l|22]. 

The standard procedure to analyze the minimal surfaces in AdS? is to consider 
the auxiliary linear problem associated with the original nonlinear equations of mo- 
tion. For any minimal surface in AdS? one can consider SU(2) Hitchin equations [6], 
which are a set of first order linear differential equations defined on the world-sheet 
with coordinates z, z. To be precise, there appear two sets of Hitchin equations 
corresponding to SL(2) L , SL(2) R evolutions of AdS?. These are promoted to a one- 
parameter family of differential equations with a spectral parameter (, where ( = 1 
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and ( = i correspond to the original SL(2) L and SL(2) R cases, respectively. The 
minimal surface is constructed as the product of the solutions of these Hitchin equa- 
tions. 

When the minimal surface possesses a 2n-gonal boundary, solutions of the Hitchin 
equations exhibit the Stokes phenomena. The world-sheet is divided, with respect to 
the asymptotic behavior of solutions, into n angular regions called the Stokes sectors. 
When moving from one sector to another at large \z\, a solution shows drastic change 
in its asymptotic behavior. This corresponds to moving from one cusp to another 
cusp at the boundary of the minimal surface. In each sector one can uniquely (up 
to normalization) choose the "small" solution, which shows the fastest decay for 
\z\ — > oo among the solutions. We let Sj(z,z;() denote such small solution in the 
j-th Stokes sector (j = 1, . . . , n). We are considering ((^-dependent) SU(2) Hitchin 
equations, and thus Sj's are 2-component column vectors. We fix the normalization 
of Sj so that 

(sj, s j+1 ) := det(s i s j+1 ) = 1. (2.1) 

Since there are only two linearly independent solutions, all s/s are expressed as 
linear combinations of any two of them. The coefficients of these linear combinations 
are called the Stokes data. The Stokes data are redundant and one finds relations 
among them. It turned out that such relations are concisely expressed in the form of 
T-system, where Stokes data are identified with T-functions [8]. In the present case, 
the T-system reads 

T+Tr = l + T^T H1 , (2.2) 

where j = 1, . . . , h — 3. The T-functions are given by 

T 2k +i(0) = (s_fc_i,s fc+ i), T 2k (0) = (s_ fc _i, s k ) + (2.3) 

for Tj with j — 0, . . . , n — 2 and the rest are set to be zero. Here we introduced the 
new variable 9 by 

C = e e (2.4) 

and the convention f^ 1 := f(9±iri/2). Note that Stokes data are by definition 
independent of z, z and depend only on (. In f!2.3p there appear Sj with j < 0, which 
are defined by successively applying the normalization condition (12.1 1) . As Sj and Sj +fl 
correspond to the same Stokes sector, they functions of z,z, proportional to 

each other 

Sj OC Sj+n- (2-5) 
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We see that T = 1 and T^_ 2 is equal to one of the proportionality coefficients 
(monodromies). The present T-system is of the standard An_ 3 type. 

While the T-functions completely characterize the shape of the minimal surface, 
they contain unphysical gauge degrees of freedom. Instead, one can consider gauge 
invariant quantities called Y-functions. In this case, they are given by 

y.i 'I J ,7,,,. (2.6) 

Note that Yj with j = 1, . . . , h— 3 are in general nontrivial functions while the rest are 
zero. Y-functions correspond directly to the physical variables. They are essentially 
the cross-ratios: 

Y 2 k = —X-k,k ,-fe-l,fc+l > ^2fe+l — — X-k-1 ,k ,_jfe_2,fc+l > ( 2 -7) 



where 



and 



._ (Si,Sj)(5fc,gl) , 
{S i ,S k ){Sj,Si) 



*WC = 1) = = : XS« , ^„(C = i) = ^ =: (2-9) 



x ik d 'jl • L ik x jl 



with The indices in the cross-ratios Xijki are labeled mod n, and the 



subscripts ± in xf, are space-time indices, not to be confused with the shift of 9. 
It follows from (I2.2p and (12.61) that the Y-functions satisfy the following Y-system 

Y+Yf = (l + Y j _ 1 )(l + Y j+1 ), (2.10) 

where j — 1, . . . , n— 3. These are the main equations that characterize the present Y- 
functions. In addition, the Y-functions for minimal surfaces in AdS% obey additional 
conditions. One is 



Y i (0)=Y j (-0) 1 (2.11) 

which follows from the reality condition of the minimal surfaces. Another important 
condition is the asymptotic behavior for small spectral parameter £ = e e — > 0, 

log Y 2k ~ ^* , log Y 2k+1 ~ ; (2.12) 

with moduli parameters Zj. This behavior is determined as follows: For large \z\, 
the small solutions behave as Sj ~ exp(w/( + where w = J z ^p(z)dz with 
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p(z) being the (n — 2)-th degree polynomial associated with the minimal surface [5]. 
By the WKB analysis one finds that the Y-functions behave as in (12.121) with Zj = 
— §^ y/p(z)dz. Namely, Zj are period integrals over cycles jj of the Riemann surface 
y 2 = p(z). 7j are taken in such a way that each of them has nonzero intersection 
with the adjacent ones j2k A 72k- 1 = 72fc A 72*1+1 = 1 [8J. For later convenience, we 
rewrite Zj as "mass" parameters 

m 2k = ~ l 2Z 2k , m 2k +i = -2Z 2k+1 /i. (2.13) 

The mass parameters are in general complex: 

m . = \m j \e* p i. (2.14) 

The other condition is the analyticity that the shifted Y-functions 

Yj(6):=Yj{6 + i Vj ) (2.15) 

are regular inside the strip 

_| <Im 0<|. (2.16) 

Incorporating these additional conditions, one can transform the Y-system rela- 
tions (I2.10p into the following TBA type integral equations 

\ogYj(6) = -\rrij\ coshtf + K jtj ^ * log(l + Y^) + K j:j+1 * log(l + Y j+l ) (2.17) 

for \<fj — <fj±i\ < vr/2, where 



1 1 

27r cosh(6 l + i(fj — i(fjt 



K SJ'W) ■= -77^77^— — T (2-18) 



and f * g '■= J f(9 — 6')g(6')d6' . The free energy associated with the TBA system is 
expressed by the solutions of the above TBA equations as 



h— 3 



" J POO JQ 

-F = A {ree = J2 — Kl cosh # log (l+F^)). (2.19) 
j =1 J -00 ^tt 

This gives the main contribution to the area of the minimal surfaces. 
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2.3. Remainder function 



The TBA equations (I2.17P completely determine the Y-functions. The area of the 
minimal surface and the reminder function are then computed by using these Y- 
functions. The remainder function is defined as 

R = A - A div - A BDS , (2.20) 

where A is the area of the minimal surface, A<n v and A-qbs are respectively the 
divergent term and the finite term read from the BDS conjecture [6]. The formulas 
for the 2n-point amplitudes are different for odd n and even n. We describe them 
separately below. 

2.3.1. Case with odd h 

In this case, the remainder function reads 

R = Ajinh + v4 peri ods + AAbds ■ (2-21) 

The first term is 

^sinh = ^frce + (H - 2)A^ , (2.22) 

where A™^ = In/ 12 is the value for the hexagon solution, which necessarily becomes 
equivalent to the regular hexagon solution. 
The second term is given by [B] 

(n-3)/2 

A vcrioAs = i (^> w ' r -^> m ' r ), (2.23) 

r=l 

where w% = $ e ^/p(z)dz and w m ' r = § mtT \J 'p(z)dz are the periods for the cycles 
with the canonical intersection form 7^; A7 m,s = 5*. In terms of the periods Zj, which 
correspond to cycles with nontrivial intersection form 9^ k = jj A 7^ (see [8]), it is 
written as 

periods = -iw jk Z j Z k . (2.24) 

Here, Wjk is the inverse of the intersection form 9^ k . 
The third term is given by [B] 

AA BDS = ^BDS-iike - ^bds = 7 £ log log ^ . (2.25) 

4 i,j=l C i,j+1 C i,j 
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^BDS-Hke is a finite term left after subtracting A div from the area of a certain reference 
region, cf^ are cross-ratios formed by nearest neighbor distances only. For example, 
when (J — i) > is odd, 

± _ X i+2,i+l X i+A,i+3 ' ' ' X j,i , 9 „n 

C i,j ~ ± ± ± ' ^z.zd; 

where Xj, Xj+i, . . . , Xj successively appear in the expression. When (j — i) > is 
even the path goes along the opposite side: Xi —> —>■•••—>■ Xj. By definition, 
c,- „■ = c,-,-. The indices of c, „■ are labeled mod n. We also define of, — 1 for \i — j\ < 1. 
For the kinematics, one can set the coordinates of the cusps so that q~. > 0. 

Note that one can express cfj in terms of Y-functions at special values by using (|2.7p . 

Case witt even n 
In this case, 

-R = ^.sinh 4" ^-periods ~\~ ^-extra + AAbds • (2.27) 

The first term A S i n h is the same as in (12.221) . The second term is [6 J 

(n-2)/2 

A pcviods = t «w m ' r -w e r w m ' r ). (2.28) 

r=2 

The third term is given by 

^extra = ~^ (w a + W B ) hgjf + ^r(w s - W s ) log 7 f , (2.29) 



where 



and 



with 



e t« s +»I> s = X 23 X 45 " " " X h,l e (w a -w s )/i = X 23 X 45 [ [ [ X h,l ^ ^ 

x 12 x 34 x n-l,n x 12 x 34 x ft-l,ra 



7l L = 7i(C = l), 7? = 7i(C = 0, (2-31) 



7i(C) = 7i(0 + 7ri). (2.32) 



Here, e ™ s+ ™ s , e^-*-)/* are given by (T^_ 2 ) {_1) " /2+1 at = -ttz/2 and 0, respectively. 
The relation between 71(C) and T\ follows from the definition of the Stokes data. We 
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see that A extra is expressed in terms of T-functions. The T-functions are expressed 
in terms of Y-functions, as we will see in section 4. 
The last term is similar to (12.25p . 



where n' = n + 1 and are defined by Cjj for the 2{n + l)-point case by the double 
soft limit x^ +1 — > xf. The explicit form is given in [22] . 

3. Integrable models and CFT perturbation 

As we saw in the last section, the problem of computing gluon scattering amplitudes 
at strong coupling is governed by TBA type integral equations. It was found [9] that 
the above TBA equations for null-polygonal minimal surfaces in AdS^ are identified 
with those of the homogeneous sine-Gordon (HSG) models [17]. In particular, when 
the resonance parameters are trivial, the HSG models admit the description of the 
bulk and boundary conformal perturbation theory (CPT). This allows us to analyt- 
ically solve the TBA equations near the CFT point in the form of high-temperature 
(or small mass) expansion. Moreover, in some special cases the TBA equations 
are also identified with those of other integrable theories. In these cases, the CFT 
analysis becomes simpler. Furthermore, from the results for the trivial resonance 
parameters one is able to analyze the case of nontrivial resonance parameters, as is 
discussed later on. In this section we review those integrable models associated with 
the TBA equations and the high-temperature expansion of the free energy. 

3.1. Homogeneous sine-Gordon model 

The HSG models are obtained as integrable perturbations |23j of the coset Gfc/U(l) rG 
CFTs (generalized parafermion CFTs) [TBI 121] ■ Here, k is the level and tq is the 
rank of a compact Lie group G. For the present purpose we focus on the case with 
G = SU(n). The classical action is 



where Sgwzw is the corresponding gauged WZW action, g is an element of G, and 



A^BDS — ^BDS-like — ^BDS — T 




(2.33) 




(3.1) 



V(g) 




tr(A + ^ f A_^) . 



(3.2) 
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A± = i\± ■ h are elements of the Cartan subalgebra f) of the Lie algebra g = su(n), 
which are parametrized by two (n — l)-dimensional vectors 

n-l 

A± = ^/i i e ±t W i . (3.3) 

i=l 

Here A, are the fundamental weights, the parameters <Tj describe the resonance when 
they are real and together with a bare overall mass scale fi, dimensionless parameters 
fli give the semi-classical masses of the solitons 



sin^ 



//* = [ijkPa, Pa = —^ (a = 1, k — 1) . (3.4) 

sm k 

The potential is identified with a linear combination of weight-zero adjoint fields. 
One can also rewrite the potential V(g) as 

2 n— 1 

v (g) = ^J2^ r M, (3.5) 

with 

/4 = Mi?'"*' » r ii(#) = ~tr((Ai ■ h)g\\j ■ h)g} , (3.6) 

so that it becomes manifest that each mass scale jUy is turned on by the field T^. 

Let us now consider the quantum theory. For the reason explained in the next 
subsection, we set Oi = 0. In the quantum theory, the adjoint fields have conformal 
dimensions 

— n 

A = A = . (3.7) 

k + n K ' 

Thus, as in the case of the sine-Gordon model |25j . on dimensional grounds the 
potential is renormalized as 

71-1 

V ~ ^2 (M 2 MjMj) 1_(A+A)/2 |"rjj(#)j . (3.8) 

As in the bare case, M is the physical mass scale and Mi are dimensionless parame- 
ters, in terms of which the physical masses are expressed as 



Mf = M lPa , Mi = MM, . (3.9) 



Replacing 



with the weight-zero adjoint operators, the action in the quan- 
turn theory may be given by 

S = kS gWZ w + A J d 2 x $ AiX , (3.10) 
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where 



= J^iMiMtf-^WikAWfo. (3.11) 

Here, <p l j are the adjoint operators which transform as the l-th (/-th) element of f) 
under the left (right) transformation. We have set <7j = and hence A = A. One 
can decompose them as 

n-l 

A = A = Mf (A+A)/2 A, , (3.12) 

where each basis vector Xj corresponds to the deformation along which Mj varies. 
Classically Xj coincide with the fundamental weights Xj, but quantum mechanically 
they are functions of the ratios of Mj. The coupling constant is related to the mass 
scale M ag§ 

A = -K n M 2 ~^ + ^ . (3.13) 

The proportionality constant K n is computed explicitly for some simple cases, as we 
will see later. 

We normalize the adjoint operators as 

(^)0r,p(O)) = « T ,|^, (3.14) 

from which 

/ \ G 2 2=i _a_ _ _a_ 

= 77^" > G(Mj) := J2 Mf +2 FijM n+2 , (3.15) 

where 

F y = A, ■ A,, (3.16) 
Note that in the classical limit Fij coincides with the inverse of the Cartan matrix. 



1 Throughout this paper, we consider the perturbation with a negative coupling constant A. 
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3.2. Correspondence with minimal surfaces in AdSs 

Let us consider the SU(n) HSG model at level k = 2. For the case with k = 2, we 
can drop the index a in Mf and absorb p a into M. In [26J, the exact S-matrix of 
HSG models associated with simply laced g is proposed. For the SU(n) 2 HSG model, 
the S-matrix reads [27] 



,s;,:,(^ = (-i)V r/t ,u,hi(0 + (7 i -<7 fe -|i) ^ (3.17) 



for j, k — 1, n— 1, where <5jfc is the Kronecker delta, Cj = ±1, and Ijk is the incidence 
matrix for su(n). Given the S-matrix, one can write down the TBA equations for 
the HSG models U7\M as 



ej (9) = MjL cosh 6 - K jk * log(l + e" £fc ), (3.18) 

where L denotes the inverse temperature and the kernel is given by 

d 

As usual, the pseudo energies €j(9) are related to the Y-functions as 

Yj = e~^. (3.20) 

It was found in [9j that the above TBA equations coincide with those for the minimal 
surfaces in AdS^ (I2.17P under the identification 

n — h — 2, MjL = \rrij\, <jj = icpj. (3-21) 

Note that the TBA equations for the minimal surfaces with complex masses 
((Pj 0) correspond to those for the HSG model with purely imaginary resonance 
parameters. Physical interpretation of such resonance parameters in the HSG model 
is not quite clear at present. In addition, when Oj — a k are nonzero, the boundary 
Yang-Baxter equations [30J are not satisfied and the boundary factorizable scattering 
may not be well-defined. These are the reasons why we have set <jj = in the last 
subsection. Also, with this restriction one can make full use of the known results 
about the high-temperature expansion, as we see in the following. We will discuss 
how to incorporate <jj in later sections. 
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3.3. Bulk perturbation of free energy 

The free energy of the model on a space of length R^> 1 with temperature 1/L (in the 
L-channel) gives the ground state energy of the model on a space of length L (in the 
i?-channel) [12J. The free energy around the CFT point is then given by evaluating 
the ground state energy of the perturbed CFT on a cylinder of circumference L with 
small coupling constant A. From the action ( 13 10p . the conformal perturbation theory 
(CPT) gives an expansion of the free energ}|f 



^free = ~C n + f n + (2tt) ^ ~ ^~ ^ 



k=l 



X 



^a,a(«i>«i) • " • ®\,\( z k,z k )) Y\(ziZi) A x dz\---dz\ 

I connected 

i=2 

(3.22) 

k=2 

Here, Z = ML denotes the scale parameter, c n is the central charge and /^ ulk is 
the bulk term. We have set z% = 1 by using translational invariance. We have 
used the fact that the coupling constant A has mass dimension 2 — 2A = an d is 

4 

proportional to I n + 2 . We have also used the fact that the one-point function vanishes. 
The central charge of the coset SU(n)2/U(l) n_1 is 

n(n — 1) . 
c„ = -L-J (3.23) 

The bulk term /^ ulk may be obtained as a generalization of the results in the 
literature [I2J|3TJ|33j|3l] • For odd n, following the procedure in [12j[33] (see also |35j) 
one arrives at the expression^ 

1 n— 1 

f^ = l* 2 T,M i (r 1 ) ij M J . (3.24) 

Note that the incidence matrix Iij has the inverse for odd n. For even n, following 
the analysis in j3l] one obtains 

/ »" ll =(^ <22 '' 2logi ' (3 - 25) 



2 We have rescaled the free energy as L 2 F/R —> F. 

3 In the course of the derivation in the references, A is assumed to be sufficiently small, so that 
the terms in the summation in (|3.22l) are smaller than I 2 . As this condition is not valid in our case 
with A = n/(n + 2), one needs appropriate modification for a rigorous derivation. We have checked 
for n — 3 that the expression (I3.24p is in good agreement with numerical results. 
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where 

n/2-l 

Q :=J2 H^+i ■ (3.26) 

3=0 

f l3.24p - fl3.26p reproduce the known results in the case of a single nonzero mass [3~T| 

SUES]. 

The coefficients /n are obtained by computing CFT correlation functions. Since 
we already know the two-point function, we find the lowest correction from CPT: 

'?S«^ ^(->-^M^). (3,7) 
Here, G is given in (I3.15P and we have used 

J d 2 z \z\ 2a \l - z\ 2b = 7r 7 (l + a) 7 (l + b)j{-l -a-b). (3.28) 
In section 6, we present the explicit expression in the case of n = 3 as an illustration. 
3.4. Single mass cases 

The TBA equations (13.181) contain n — 1 mass parameters. When some of them 
are set to be zero, the corresponding massless pseudo energies 6j(9) become finite 
constants for L — > oo, while massive ones become infinite [3 1 J . This gives rise to 
the appearance of CFT with a different central charge. Nevertheless, one finds little 
difference concerning the small mass expansion: When one turns on/off some of the 
masses, only the constant term changes discontinuously and the other expansion 
coefficients changes continuously. This can be checked numerically. 

For HSG models with general masses, the precise form of the perturbing operator 
(13.111) . as a function of masses, has not yet been well understood. On the other hand, 
when only one mass parameter is nonzero, the CFT in the small mass limit becomes 
simple and the perturbing operator as well as the exact coupling-mass ratio is known. 
This allows us to explicitly compute the expansion coefficients at low orders. Below 
we present a list of such cases and collect useful facts in performing the small mass 
expansion]^] Remarkably, combining the cases in subsections 3.4.1 and 3.4.3, we are 
able to analyze a case of two general masses, which we will discuss in section 6. 

4 The reader is referred to 37 for a more extensive classification of conformal perturbations 
associated with TBA equations. 
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3.4-1 ■ Integrable perturbation of unitary minimal model 
Let us first consider the case 

Mi = M, others = 0. 



(3.29) 



The TBA equations (13.181) in this case are identified with those of the (RSOS)„ 
scattering theory [3T1[32]. The (RSOS)„ scattering theory is regarded as the massive 
perturbation of unitary minimal model M. n +i,n+2 by the primary field $13 with 



dimension A = A 



n+2 ' 



The action for the (RSOS) n scattering theory is 



S(RSOS) n 

where the relevant operator is normalized as 

$1,3(^1,3(0) 



SMn+i,n+2 + ArsOS / dfx $1,3 



1 



I4A 



and the coupling Arsos is related to the mass Mi as 

Arsos 

with 



. K RSOS M 2(l-A) 



RSOS 



1 (n + 2) 



7 



3(n + l)\ fn + l 



n + 2 



7 



n + 2 



v^r(^; 
2r(2±i) 



4 

n + 2 



nn(2n + 1) 

and j(x) = r(x)/r(l — x). 

Taking into account the normalizations A3. 15j) . (13.31 j) . one finds that 



(3.30) 



(3.31) 



(3.32) 



(3.33) 



XG(Mi 



Af<=0(#l) 



A 



RSOS 



and hence 



A = -k„M™+2. 



-RSOS 



(3.34) 



(3.35) 



3-4-2. Integrable perturbation of unitary SU(2) diagonal coset 
The above case is generalized to the cases 

M fc = M, others = 0, for jfe = 1, . . . , n - 1. 



(3.36) 



It is known that the TBA equations (I3.18P in these cases describe the system obtained 
as integrable perturbation of (SU(2) fc x SU(2) n _ fc )/SU(2) n coset CFT [3SIEZ]. The 
perturbing operator is the primary field 0i,i,3, which corresponds to the branching 
of the product of two trivial representations into the adjoint representation. This 
operator has conformal dimension A = A = The exact coupling-mass ratio in 
these cases is found in |38|. 
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3.4-3. Integrable perturbation of non-unitary SU(2) diagonal coset 

Let us next consider the cases where 

M k = M n _ fc = M, others = 0, (Jfe = 1, . . . , n - 1) (3.37) 

with n being odd. This kind of configuration is invariant under the Z 2 outer automor- 
phism of the A n _i Dynkin diagram. In this case, one can regard the TBA equations 
as those corresponding to the tadpole diagram T( n _i)/ 2 (= A n _i/Z 2 ) with a single 
mass parameter being turned on. It is known that they describe the system obtained 
as integrable perturbation of non-unitary (SU(2) fc x SU(2) n / 2 _ fc _ 1 )/SU(2) n / 2 _ 1 coset 
CFT [37] . The perturbing operator is 01,1,3 with dimension A = A = 
In particular, the case with k = 1, namely 

Mi = M„_i = M, others = (3.38) 

corresponds to integrable perturbation of the non-unitary minimal models M. n ,n+2- 
The case with n = 3 and k — 1 will be used in the analysis for the decagon in 
section 6. 

4. Conformal perturbation of g- and T-functions 

As we saw in the previous section, the relation between the free energy in the L- 
channel and the ground state energy in the i?-channel allows one to derive an ex- 
pansion of the free energy near the CFT limit. Such an expansion is studied for 
the 6-point amplitudes with the AdS§ kinematics in |13j . To obtain the full ex- 
pression of the amplitudes, one further needs the expansion of the cross-ratios or 
the Y-/T-functions. A key observation [39| l4*0] for this purpose is the relationship 
between the T-functions and the ^-functions |14j . In this section, we discuss the con- 
formal perturbation of the T-functions of the HSG model associated with the coset 
SU(n) 2 /U(l) n_1 . We follow [1~5||4"0] where the perturbation of the g- and T-functions 
in the ADET purely elastic scattering theories is discussed. The discussion below 
holds for general n with trivial formal monodromy (i.e., T n = 1) and will be applied 
to concrete examples in the following sections^] 

4.1. (/-functions in homogeneous sine-Gordon model 

We start by considering the partition function Zi a \ a \ [L, R] on a cylinder of circum- 
ference L, length R, and boundary conditions of type a at both ends. It is expanded 

5 For even n, one can consider the cases with non-trivial T n , where the discussion may need some 
modifications. Such cases appear for 2n + 4 > 12, which we hope to discuss elsewhere. 
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by the eigenvalues of the circle Hamiltonian H cac (M, L) as 

oo „ 

Z {ala) [L,R] = (a\ e-*^W)| a) = e'™?^ , (4.1) 

where M is the mass scale defined through = 



with \ip p ) being the eigenstates of the Hamiltonian, and I = ML. The g-function is 

ii/i/l 1 \t - iiKli'.w.fiiwv -\- 1 i / \ 1 i i i / \ • i t - 4 /n - i i i ill / i i -, \ i i i ( , , 



then defined by subtracting the linear term in L from Gw,(l)' 



\ogg la) (l) = \oggl° Q ) ) (l) + f lQ) L. (4.3) 

f\ a ) is the constant boundary contribution to the ground state energy Eq V1P of the L- 
channel Hamiltonian H^£JR). The (^-function counts the ground state degeneracy 
and is known to decrease along the renormalization flow [14J. It is also a subleading 
contribution to the partition function for large R. 

To consider the ^-functions in the HSG model, we assume that the HSG model 
admits an integrable generalization with boundaries, and the boundary scattering 
amplitudes are also diagonal. In this section, we also set the parity breaking param- 
eters <Tj = as discussed in section 3. We will discuss how to incorporate Oj in later 
sections. 

In the presence of boundaries, one has the boundary reflection factors Rj(0), 
which are constrained by unitarity and crossing-unitarity [4"T||4"2]. 



Rj(9)Rj(—8) = 1 , Rj(9)Rj(e - iix) = 5^(20) , (4.4) 

where the anti-particles are the same as the particles in our case, J = j. The boundary 
Yang-Baxter equations are indeed satisfied for ajk :— <7j — = 0. In general, 
reflection factors also have to satisfy the boundary bootstrap equations. They are, 
however, trivial in our case, since the bulk S-matrix (I3.17P with <jjk = does not have 
poles in the physical strip < lm.6 < it. The form of the constraints (I4.4p shows that 
a set of reflection factors Rj(0) generates another R'j = Rj/Zj, if Zj(8) are a solution 
to one-index versions of the bulk unitary and crossing-unitary equations 



Zj(9)Zj(-e) = 1 , Zj(9)Zj{0 -m) = l. (4.5) 
In particular, we use a solution, 

Zf c \9):= ((1 + CV1 " cVf , (4.6) 
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in the following, where 



sinh M9 — i^x) 



A boundary a is then associated with a set of the reflection factors Rj (9). It 
turns out that the g-f unction corresponding to |a) satisfies the following integral 
equation PSIHSHHSS], 

log g\ a) (l) = logC| a) + S(0 (4.8) 

n-1 



it,— ± n. 

+ 69 (0i Q> W-^)-20 ii (20))l O g(l + e^W) 



Here, C\ a ) is a symmetry factor associated with the vacuum degeneracy at infinite I. 



'' and (pjk are given by the boundary and bulk S-matrices as 



cj)f{9) = -de\ogR\ a \9) , <f> jh (9) = ^d logS jk (9) . (4.9) 

7X1 Z711 

S(Z) is a certain boundary-condition independent term, precise form of which is 
irrelevant for our purpose. For details, see [T5t fT6]. 



4.2. Relation between g- and T-functions 

In the conformal limit, boundary conditions are labeled by primary fields, and hence 
a boundary a may be specified by the corresponding primary field. To describe the 
relation between the g- and the T-functions, we then consider the reflection factors 
corresponding to the boundary condition labeled by the identity operator, R J 1 \9), 
together with the deformed ones, 

<' C >W = (4.io) 

R^ is expected to be minimal, namely, having the smallest number of poles and 
zeros. We assume the existence of R^P • These reflection factors give the ^-functions 
gu\ and g\k,c) through f)4.8p . the ratio of which satisfies an integral equation, 

log ( (4.11) 



9\i)(l) C|fe,c), 

/ 47 (cosh(0 + ifC) + cosh(0-zfC)) l0g ( X + Yk{ -°\ 
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We notice that almost only the information of the deforming factors Z^' ^ has re- 
mained. 

On the other hand, from the T-system (12.21) and the relation between the T- and 
Y-functions (I2.6p . the T-functions obey the integral equation, § 

\ogT k {6) = -i/ fc cosh0 + tf*log(l + y fc ). (4.12) 

Here, the kernel is given by (I2.18P with (pj = ipji = 0: K = l/(27rcosh6 l ). v k 
specify the asymptotic behavior of T k {6) and are related to the mass term of the 
TBA equations as 

m k = u k -i + Vk+i ■ (4.13) 
Comparing (14.111) and (I4.12p . and using T k (0) = T k (—9), one finds that 

log ( = ^cos&U logli^C) . (4.14) 



9\i)(l) C\k,c)J \2 J V 2 

Since C\ a \ is associated with the vacuum degeneracy, which may be determined by 
the symmetry, we expect Cm = C\ k; c)- Assuming this and subtracting the linear 
term in / oc u k , we arrive at an important formula, 

C (0) 



TAi-C). (4.15) 



\k,C) 

Mo) -«V"2' 
4.3. Expansion of T-functions 

Using (I4.15P and the conformal perturbation, one can derive an expansion of T k {9) 
near the CFT limit. To this end, we first note that the conformal perturbation gives 
the expansion of the ^-functions [EHUD], 



log0, ( g(A,L) = X>l 1} ( AL 



2-2A\<j 



9=0 

oo 

logSg(A,M,L) = ^^(^- A ) P (AL 2 - 2A )^, (4.16) 

p,g=0 

where /i is the coupling of the boundary perturbation, and we have assumed that 
the dimension of the boundary perturbing operator is A = njin + 2). Following the 



For even n, an appropriate gauge has to be chosen. 
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argument in [T5], we also assume the relation between the boundary coupling and 
the deformation parameter, 

fi = fi cos( — — C) , (4.17) 

where /j,q is some constant. Together with (|4.15|) . this relates the boundary coupling 
/i to the argument of T^{9). (14.1 7p also means that the boundary coupling vanishes 
at C = (n + 2)/2. In this case, the boundary conditions become conformal and are 
described by (linear combinations of) the Cardy boundary states. The conformal 
perturbation with only the bulk coupling A turned on then gives [T5|l4T)] 

log£| Q) (A,L) := log(a|fi) = hgg {a) + Xd[ a) L 2 ^ + ■■■. (4.18) 

Here, \fl) is the full ground state, 

and B(a,b) = T(a)T(b)/T(a + b) is the Euler beta function. We have also introduced 

g la) := (a\0) , gf a) := (a|$) , (4.20) 

with |0) and |$) being the CFT vacuum state and the state corresponding to the 
bulk perturbing field $ = $ A x, respectively. Q\ a ) is different from in (14. 2 p by a 
normalization factor (filfi) 1 / 2 , but this factor is cancelled in the formula (I4.15p . 

On the T-function side, we first note that the Y-functions for the 2(n + 2)-point 
amplitudes have the periodicity, 

/ n + 2\ 

Yk(e + in——) = Y n _ k (6). (4.21) 



2 

Yfc are also analytic for finite 9, and even functions of 9, Yk(9) = Yk{—9). These 
properties are common to the T-functions, which leads to the expansion 

oo _ 

Tk{9) = E c l P) (0cosh(-^^), (4.22) 

p=0 

with 

c? q) (l) = c™(l), cf +1) (0 = -et 1} (0 (geZ> ). (4.23) 

For small I, the coefficients behave as c^f (I) ~ /( 1_A ^ P , since the Y- and T-functions 

show plateaus for — log(l/7) <C 9 <C log(l//). The conformal perturbation (14.161) 
further suggests that cf\l) are expanded as 

oo 

c ( *\l) = J2$' 2q) l (1 ~ A)(p+2q) - (4-24) 

9=0 
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(See Appendix C and D for details.) The Y-functions have a similar expansion. 
Similar double expansions have been discussed also for other TBA systems [13,46]. 
Substituting this expansion into the Y-system (I2.10p . one then finds that some lower 
coefficients vanish. In term of the T-functions, the result reads 

4 0) (0 = 4 ' 0) + o(/ 3 ^)), ci 1) (o = o + o(/ 3 ^)), 

c f(l) = c f> °)/2(l-A) + ( /3 (l-A)) (4 25) 

Now, from (14151) . f HlSj) . fl422|) and (1Q5]) with C = (n + 2)/2, the expansion of 
Tfc(#) is determined. Comparing both sides of (I4.15p . one first notices that the two 
expansions (14.161) and (14.251) are consistent with each other, once = is taken into 



account, meaning that the one-point functions of the boundary perturbing operator 
vanish in unitary theories. Furthermore, 

T k {9) = 4°' 0) + c j&°) Z 2(i-A) cosh(-i^) + 0(l^) , (4.26) 



where 



(0,0) _ 9\k) (2,0) _ 9\k) ( ,\k) 



°k ~ „ ' °k 



^(<>-<> n , (4.27) 



9\t) ' 9\i) 

and we have used (I3.13p . We have also set g\k,c) ='■ 9\k)i fl 1 * c) =: ^fl) anc ^ =: di \ 
since they are evaluated for the unperturbed boundary states which are independent 
of fi and hence of C. 

In the course of deriving the above formula, we have made several assumptions 
following [151 HQ]: the existence of integrable boundary perturbations of the HSG 
model by operators with dimension A, that of the reflection factors associated with 
the identity operator, the invariance of the symmetry factor C\t) under the deforma- 
tion, and the relation (I4.17p . We will check that these are consistent with numerical 
computations and the results in the conformal limit, which we discuss shortly. 

4.4. Identification of boundary conditions 

When evaluating the expansion (I4.26p . we need to identify the boundary conditions 
represented by the reflection factors (I4.10p . For this purpose, we first recall that the 
Cardy boundary states are of the form 

l a > = E4f=H' (4-28) 

o V^Op 
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where \p)) are the Ishibashi states. In our case of the SU(n)2/U(l) n_1 coset CFT, 
the modular S-matrix, S pp /, is given by the product of the S-matrix for SU(n)2 
and the complex conjugate of the S-matrix for U(l) n_1 : S pp > = S^S*^. Here, 
p := [pi, p n -i] is the Dynkin label of su(n). Since we deal only with primaries 
whose u(l) weights are zero, the u(l) part is trivial and hence we drop S*^ in the 
following. 

In general, the index of the T-functions labels the representations of the un- 

Ifc c) 

derlying symmetry. Thus, we infer that Rj ' correspond to the k-ih fundamental 
representation with the Dynkin label whose components are (pjjj = ftj.k- It then 
follows that 



9\i) 



S 



S, 



(2) 



9\k) 



(2) 



S, 



(2) 



(4.29) 



^oo y ^oo 

Here, the S-matrix for SU(n) 2 is given by the formula 



C(2) 



(n + 2) 



-(n-l)/2 



jn{n-l)/2 



n 



x det exp 



2m 



exp 



'n-l 



n + 2 




3 (to + !) 



(4.30) 



l<a,b<n 



From fl4T29]) and (tODjl . we find the CFT limit of the T- and Y-functions: 

(fc + l)7T 



(2) c (2) 



c(2) 

(2) 
00 



sin 



n+2 



sin 



n+2 

&71 



sin sm 

n+2 



(fc + 2)7T 

n+2 



(^00 ) n+2 

These agree with the result in [8], which supports our formula (I4.26P and identifica- 
tion of the boundary conditions. We note that the relation between the (/-functions 
and the T-functions at the CFT point naturally explains the fact that the quan- 
tum dimensions (ratios of the modular S-matrices) are the solutions of the constant 
T-system, called the Q-system [TUISSISS] - From the relation between the T- and 
(/-functions, this may hold for general TBA systems. 

Since the bulk perturbing operator $ is a linear combination of the adjoint oper- 
ators, are similarly expressed by the elements of the modular S-matrix with the 
Dynkin label of the adjoint representation p ad y 



<i> 
9\i) 



G(M % 



o(2) 
°0p adj 



c(2) 
°0p adj 



<l> 

9\k) 



G{Mi) 



o. 



(2) 

PfcPadj 



c(2) 



(4.32) 
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Thus, at the leading order, both of the expansions of the free energy and the Y- 
functions are given in terms of G(Mj). In addition, in order to obtain the full 
expressions of the scattering amplitudes, we need the explicit form of $, as well as 
the bulk coupling A, in terms of the TBA masses m&. These are discussed in section 
6. 

5. Remainder function for the octagon 

As we have seen in section 2, the minimal surface with a null polygonal boundary in 
the AdS space is described by the Y-system or the TBA equations. The Y-functions 
and the T-functions play important roles in this picture. In the previous two sections, 
we have seen the relation between the (^-functions and the T-functions in the under- 
lying integrable model. This relation enables us to compute the high-temperature 
(small mass) expansions of the T-functions by using the conformal perturbation 
technique. Consequently, the remainder function is expanded around the kinematic 
configurations associated with regular polygons. Here we consider the first non-trivial 
example in AdS^. the octagon. In this case, the exact expression of the remainder 
function at strong coupling has been computed by Alday and Maldacena [6], and one 
can learn much about the expansions of the T-functions and the ^-functions. The 
underlying integrable model corresponding to the octagon is the off-critical Ising 
model (with a complex mass). The exact (^-function for the off-critical Ising model 
was obtained in P^ | IqT7 | [HT] . Though the TBA system of the off-critical Ising model 
is trivial, the high-temperature expansions of the free energy and of the T-function 
are still non-trivial. 

As in (I2.27p . the remainder function is divided into several pieces: 

7tt 

-^8 = "7T + ^frcc + ^periods + ^extra + Av4 B DS- (5.1) 


These terms are given by 

r°° a+ 

A iree = / — Zcoshtlog (1 + e - lcosht ) , (5.2) 

^■periods = 0, (5-3) 
^extra = - ^ (COS (f) log 7^ + SU1 (f) log ) , (5.4) 

AAbds = ^BDS-iikc - A BDS = ~ log(l + X ~) log (l + ^J, (5.5) 



where 



X+ = e lsiD t X - = e~ lcos t (5.6) 
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and 7^,7i are obtained from (I2.3ip with 

log 7 i(C = e e ) = -L , log (1 + e-' cosht ) . (5.7) 

Z7T J_ 00 cosh(t — + 1(f)) 

We note that (I5.7P is a special case of (I4.12p for k = 1 with z/i = and Yi(9) = 
e _ ' cosh£l , in accord with (12.321) . Our goal here is to expand the remainder function 
around / = 0. For the octagon, we can get the all-order expansion at arbitrary 0. 

The expansion of the free energy was studied in [33] , and the result is the follow- 
ing, 

A free = ^-^(logy + ^ + log 7 r-7E 

+-E^iiH i -^iV(2^ + i)^y fe+2 , (5. 



k=l 



where 7e is the Euler constant. 

In order to expand A extra , we introduce the following function, 

where 7]" and 7^ are related to this function as 

log 7l L = F(/,0), log 7l R = F (l,<l>- |) . (5.10) 

Note that this function is related to the exact (^-function for the off-critical Ising 
model considered in [16]. In [16], the exact (^-function was discussed in two special 
cases, which essentially reduce to ip = in our case, ip 7^ corresponds to the case 
that the boundary magnetic field is turned on in the off-critical Ising model. Here we 
obtain the small / expansion of the g- function for general values of ip with \(p\ < 7r/2. 
As we will discuss in Appendix A, this function has the expansion 

F(1,V) 

= il„g(l + e-<- )-JL 



cos <p> ^log j + 1 + log 7r - 7 E ^j + (p> - ~) sin ip 

+ £ (?) sVt (' - 2^) a2k + r) cos " ^ { - k ' 1; I 8in2 *> © 



2fc+l 



(5.11) 
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Using this expansion, we find the small I expansion of A extra , 
A extra = -^cos01og(l + e- M ) - ^sin01og(l + e~ Uo ^) 

I 2 f I , 7T 

+ — log - + 1 + log 7r - 7 E - - cos <p sin . 

47T \ I 1 

-tidM 1 -^) + y^AW f^f" (5.12) 



fc=l v / \ / 

where the function fk{4>) is expressed in terms of the hypergeometric functions, 
fk{4>) = cos 2 2 Fi(-k, 1; - - k] sin 2 0) + sin 2 2 Fi(-k, 1; ^ - /c; cos 2 0). (5.13) 

Combining ( 15. 5p . ( 15. 8 p and (15.12p . we obtain the expression, 

57T 1, / /cos0\ / , lsm6\ I 2 
R S = - - - log ^2 cosh — j log ^2 cosh — j + - 

i \ , , / k + i rt \ rr 2k+2 



7T 



(5.14) 



Note that the non- analytic terms in ( 15. 8 p and ( I5.12p cancel each other out. The 
expression ( I5.14p is convergent for \l\ < it. One can immediately confirm that the 
remainder function is expanded in I 2 : 

oo 

Rs = J2 R s k) ^ l2k > ( 5 - 15 ) 

fc=0 

where the first four coefficients are given by 

(0) 5vr log 2 2 
#8 W = "4 > ( 5 - 16 ) 

^(^^ + (^r-^)co S 4, (, 18 ) 

, 6) = -81og2 + 3 _ f 81og2 + 5 _ jlggA 

8 73728 V 122880 12807T 5 ; v y ' 

We note that these maintain the symmetries </> — > —(f) and <p — > <fi + |, which are 
due to the space-time parity and cyclicity [B]. The expansion ( I5.12p was derived for 
< < 7r/2, but is valid for arbitrary due to these symmetries. In section 7, we 
will compare these results with the remainder function at two loops. 
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6. High-temperature expansion for the decagon 

In the previous section, we have considered the remainder function for the octagon, 
and have obtained its all-order expansion with respect to the mass scale parameter 
I. The crucial point there is that the TBA system for the octagon is trivial. For the 
general 2fi-gon (n > 5), however, analytic solutions of the TBA equations have not 
been known yet. In this section, we consider the second simplest case n = 5: the 
decagon, and see how to compute the high-temperature expansion of its remainder 
function. Here, the underlying integrable theory is the homogeneous sine-Gordon 
model associated with the SU(3) 2 /U(1) 2 coset CFT. 

The remainder function for the decagon is divided into the following parts, 

7 

-RlO = + ^periods + ^frec + AA B DS- (6.1) 

The period part is given by 

^periods = —g{miffi2 + m 2 rh 1 ) = -\^M x M 2 l 2 cos(y?i - (p 2 ), (6.2) 
where rrtj = MjLe tlf,J , Mj = Mj/M and I = ML. The free energy part is written as 

A iiec = Y, / —MjLcosheiogil + Yjie)), (6.3) 

where Yj(9) is defined in f)2.15p . From f )2.25p . the part Av4 B ds is given by 

AAbds = \ log -$- log (6-4) 

i,j=l C i,j+l C i,j 

For the decagon, there are four independent cross-ratios. This is consistent with the 
fact that the TBA system has two independent complex parameters rrij (j = 1,2). 
The cross-ratios cf 3 and cf 4 are related to the Y-functions, 

(7TZ \ 
-yj , (6-5) 

c+ = F 2 (0), q 4 = Y 2 . (6.6) 

The other cross-ratios cf 4 , cf 5 and cf 5 are expressed by and cf 4 , 

± _ 1 + 4 ± _ i + 4 ± a . i \ A . i 

c 24 — + ' c 35 — + i c 25 



C 13 C 14 V C 13/ V c 14 



1 + — 1 + — " !■ (6-7) 



Our goal is to find the small / behaviors of the remainder function. 
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6.1. Case with real masses 



Let us start by restricting our attention to the case that two masses are real: ipi = 
(p2 = 0. In this case, we can directly use the results of the CPT in section 4. The 
period term is 

^periods = ~MlM 2 l\ (6.8) 

Let us consider the free energy part. Since the central charge of the SU(3)2/U(1) 2 
coset CFT is C3 = 6/5, the free energy goes to 7r/5 in the limit (1711,1712) — >■ (0,0). 
The perturbing operator has the dimension A = A = 3/5. From (I3.24p . the bulk 
term is given by 

/ 3 bulk = X -M x M2l 2 . (6.9) 
Thus the free energy has the following expansion, 

A iree = ^ + \m x M 2 1 2 + £ fi k) l ikl \ (6.10) 

k=2 

where /g is computed by the CPT (see (13.221) ). In particular, is read from 
dS2ZD: 

f? = l^G\M u M 2 )Cf\ Cf = 3(2 7 r)t 7 2 g) 7 (-i), (6 . n) 
where j(x) = T(x)/r(l — x) and 

2 

G(Mi, M 2 ) = M^FijMf 5 . (6.12) 

We need to determine the symmetric matrix F, which has two independent compo- 
nents for the decagon. In principle, it should be possible to fix them by considering 
the quantum theory of the SU(3) 2 /U(1) 2 HSG model. However, we take a different 
strategy here. Fortunately, we can completely fix by the following consideration 
in the case of the decagon. 

Let us first take the limit (Mi,M 2 ) — > (M, 0). In this limit, the TBA equations 
for the decagon reduce to those for the (RSOS)3 scattering theory. Therefore the 
correction (16. lip should be equal to that for the (RSOS) 3 scattering theoryjl and we 



7 Since the (RSOS)3 model has the central charge c = 7/10, the free energy goes to 7n/60 in 
the limit ML — > 0. The difference of the constant terms of the free energies in two theories comes 
from the fact that the contribution of particle 2 is absent in the (RSOS)3 scattering theory. One 
can numerically check that in the homogeneous sine-Gordon model, the contributions of particles 
1 and 2 go to 77r/60 and 7r/12, respectively in the limit ML — > with M2 <C Mi = 1, and the 
sum of two-particle contributions gives the correct value 7r/5. Furthermore all the high-temperature 
corrections of particle 2 should vanish if we take the limit M2 0. 
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obtain 



F u k 3 = 4 sos , (6.13) 

where nf sos is given by (I3.33p . 

Let us next consider the case: Mi = M 2 = M. In this case, the TBA equations 
are regarded as those associated with the tadpole Dynkin diagram T\. The corre- 
sponding integrable model is the non-unitary (SU(2)_ 1 / 2 x SU(2) 1 )/SU(2) 1 / 2 coset 
model perturbed by the primary field 0i,i,3 with dimension A = A = 3/5 (see sub- 
section 3.2). We discuss the perturbation of the above coset model in Appendix B. 
Using the result (IB. 101) . we find the constraint 

^ 4(Fll + Fl2M ^U^%(T,(-l)^(%(t)- (6-14) 



3 iK 11 1 1Z > A 8V47 'Vv ' v sy ' vsy ' Vs 

Combining (16 . 1 3[) and (16. 14p . we obtain 

F a 1/3 V /5 /l^ 1/s 



1+ tr-2{v>) HsJ • (615) 

From these two considerations, the order Z 8 / 5 correction of the free energy must have 
the following form 

,(2) _ ,(2) p-2^2 



where 



and 



fr = f^os ) Fu'G'(M l ,M 2 ) (6.16) 



F^-G(Mi, M 2 ) = M x 4/5 + M 2 4/5 - BM^M^ 5 , (6.18) 

As seen in section 3, the matrix F deviates from the inverse of the Cartan matrix. 
We have confirmed that (16.161) is in good agreement with the numerical results for 
arbitrary M± and M 2 (Fig. [Q. 

Now we proceed to the expansion of AAbds- in order to know the small I behavior 
of AAbds, we need the high-temperature expansions of the Y-functions. Since the 
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Figure 1: (a) The scale parameter /-dependence of the free energy with fixed 
M1/M2 = 1, 1/2 and 1/10. Dashed lines represent the curve f + \M x M 2 l 2 + / 3 (2) / 8/5 . 
Deviation from the analytic formula at I = 0.5 for Mi/M 2 = 1/10 comes from 
the next order correction 0(7 12//5 ), which is estimated from the numerical fit as 
0.08Z 12 / 5 ~ 0.015. (b) Mi-dependence of the coefficient of l 8 ^ 5 of the free energy 
for M 2 = 1- Dashed line corresponds to the curve given in (16.161) . 



Y-functions are related to the T-functions, we can use the results in section 4. In 
the decagon case, the relations between Yj and Tj are as follows, 



Y 1 {6)=T 2 {6), Y 2 {6) = T l {6). 



(6.20) 



As in (4.19) and (4.27), the high-temperature expansions of the T-functions are 
computed from the data of the g- functions at the CFT point. Using the formula 
(4.30) of the modular S-matrix, one obtains 



M = M = 2 cos(^V 9 ^ = G(M 1 ,M 2 ) 
9\i) 9\i) v 5/ g w 



71 

2 cos I — 

5 



1/2 



^>=^t = -G(M 1 ,M 2 ) 
9\i) 9\2) 



71 

2 cos ( — 



-3/2 



Therefore from (I4.26p . we find 

Yj(9) = + y< 2 >(M ls M 2 )Z 4 / 5 cosh (j\ + o(/ 6/5 ), 



where 



Y (0) = 2 cos ( 



(6.21) 
(6.22) 



(6.23) 



(6.24) 
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and 

.(2) 



Y^{M U M 2 ) = yf^ os) F^G{M x ,M 2 ) (6.25) 

«sos )3 - ^ ('I) [lO «» (f ) 7 (|) 7 g) ] V2 • (6.26) 
From (16. 5p and (16. 6p . the cross-ratios are expanded as 

c - = c + = y(o) + y(2)/4/ 5 + 0(/6 /5 ); (627) 



-13 — °14 



As shown in Appendix D, the expansion of the AAbds is largely constrained by the 
Y-system, the structure of the conformal perturbation, as well as the symmetries 
associated with the space-time parity and cyclicity, under which AAbds is invariant. 
Consequently, it turns out that the terms of C(/ 4//5 ) are enough to give the expansion 
of AA B ds U P to (D(l 12 / 5 ). We then obtain the high-temperature expansion of AA B dSj 

AAbds = -^log 2 (2cos (I)) + B 2 (Y^fl^ + 0(/ 12 / 5 ), (6.29) 
where B 2 is given by 

B2 = 20 cos< log (2cos (f))). (6.30) 

Note that the first order term (9(/ 4/ ' 5 ) vanishes. In summary, from ( 16. 8p . (I6.10p . (I6.16P 
and (I6.29p . the remainder function with the real masses has the following expansion, 

i? 10 = <+ J RS 4 ) Z 8/5 + O(/ 12/5 ), (6.31) 

where 

(6.32) 



(0) 


39 

7T - 

20 


10 — 


(4) 


(4 


10 — 





2 V V5 



4? = (-^tan (|) + B 2 ) F( 2 )(M 1; M 2 ) 2 . (6.33) 



Note that A per i 0( j s is canceled by the bulk term in the free energy, and the remainder 

f(2) // (2) 

'(RSOS) 3 / ^(RSOS) 3 ' 



function is expanded in / 2//5 . We also comment that the ratio /(rsos^^rsosyJ 2 
interestingly becomes very simple, 



f (2) 1 

■/(RSOS) 3 1, fK\ , aQ A\ 

— ^ = -5 tan UJ- (6 - 34) 

^(rsosv 
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6.2. Case with complex masses 

So far, we have focused on the case that two masses are real. We would now like 
to discuss the general situation where two masses are complex. The phase of the 
complex mass corresponds to the purely imaginary resonance parameter. As dis- 
cussed in section 3, it is not clear if this case can be treated within the framework 
of the conformal perturbation of the HSG model. However, one can expect that the 
expansion is analytic in the mass parameters [16], and the expansion for the real 
masses can be extended to that for the complex masses by continuing the mass pa- 
rameters. This is also expected from the point of view of the TBA equations. In fact, 
we will see that the results obtained in this way are in agreement with numerical 
computations. Furthermore, one can arrive at the same conclusion for some relevant 
quantities by considering the chiral limit of the TBA system, which is discussed in 
detail in Appendix C. 

The way to incorporate the phase is determined by the the following facts: the 
resonance parameters in the TBA equations are understood as due to the rescaling 
of the mass parameters, the free energy, by definition, should depend only on the 
difference of the phases ipu = <f \ — <f2, and A, A are of the form (13.31) semi-classically. 
We thus make a replacement, in the complex mass case, 

A -> YsiMje^y-^+^Xj, (6.35) 

3 

A -> ^(M ie -^) MA+A)/2 Ar (6-36) 

3 

Then the two-point function of the perturbing operator becomes 

(* x,x(z)$ a,a(0)> = 1 N ; 2/5 (6.37) 

where G is given by (I6.18p . 

Let us consider the free energy. As mentioned above, the free energy must be a 
function of ipi 2 - This suggests that the bulk term is modified as 

1 2 1 

/3 Ulk -> / ^(Mie^ir 1 )^^) = -M.m 2 cos cp 12 . (6.38) 

Taking into account these modifications, we find that the expansion of the free energy 
is given by 

A bee = ^ + ^M X M 2 / 2 cos V12 + f^sos^nHM^^^e^l 8 / 5 + 0(l 12 ' 5 ). 

(6.39) 
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Figure 2: (a) The /-dependence of the Y-function Yi(0) for fixed Mi/M 2 = 1, 1/2 
and 1/10 at (fx = (p 2 = n/20. Dashed lines represent the curve (16.401) at 9 = 
up to the order Z 4 / 5 . Deviation from the analytic formula comes from the next 
order correction 0(Z 6 / 5 ), which can be estimated from the numerical fit as 0.15/ 6 / 5 ~ 
0.004 (/ = 0.05) for M x = M 2 = 1. (b) Plots of the coefficient of l A ' b in Ki(0) for 
ipi = 7r/20 and various tp 2 a t 2Mi = M 2 = 1. Dashed line represents the curve 
i (F( 2 )(M ie -^,M 2 e-^) + Y^(M ie ^,M 2 e^)) =: y$ sosh h(M v ^). 

The relation between the g- and T-functions in section 4 is not applied to the 
case of complex masses. However, by similarly complexifying the mass parameters, 
we obtain the expansion of the Y-functions, 

1 



Yj(0) = 2 cos 



7T 



M 2 e-^ 2 )e 4e/5 



+ 



y( 2 )(M ie - 

Y^\M ie ^\ M 2 e^ 2 )e- 49/5 V /5 + £>(Z 6/5 ) 



(6.40) 



where is given by (16.251) . We have checked that this formula agrees with the 
numerical results (Fig. |2]). The expansion of the Y-functions for the complex masses 
is also discussed in Appendix C from the chiral limit of the TBA system. The 
space-time cross-ratios are again obtained by using the relations (16.51) and (16. 6p . In 
addition, using (I6.40p . the expansion of AAbds is given as in the case of the real 
masses by 

'7P 



AA 



BDS 



5 1 2 

— log 

2 5 



2 cos 



B 2 \Y {2 \M ie ^\ M 2 e^ 2 )\ l 8/5 + 0(l 12/5 ). (6.41) 



Collecting all the above results, we finally find that the remainder function with 
the complex masses has the expansion, 



i?io = M°o ) +i?S 4 ) / 8/5 + 0(/ 12/5 ), 



(6.42) 
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Figure 3: The /-dependence of the remainder function with (a) equal phase (p± = 
ip 2 = tt/20 (b) different phase ip± = n/20, <p 2 = 7r/5. Dashed lines represent the 
curve R^+R^l 8 / 5 . 



>(0) 39 5 , o /„ /7T 



where 

Rio = ("^ tan (|) + B *) |^ (2) (Mie^,M 2 e^ 2 )| 2 . (6.44) 

F (2) and B 2 are defined by (16.25P and (I6.30p . respectively. Note that the A per iods is 
canceled by the bulk term of the free energy again. We have confirmed that this 
formula for R\q is in good agreement with the numerical results for various values of 
Mi, 2 and (p lj2 (Fig. [3]) . 

At the end of this section, we comment on the relation between cross-ratios and 
the parameters in the TBA system. In order to express the remainder function as a 
function of the cross-ratios, one has to invert the relations (I6.5P and (I6.6p . This is 
complicated for general complex masses. However, when the phases of rrij are equal, 
i.e., ipj = tp, Yj(9) are obtained from those for the real masses by the shift 9 — ^ 9 — ip. 
One then simply has 

U = tan- ( cot ( Vi P^f] , YW* = ^ , (6.45) 

5^ V V 5 / S^ + Sc^J cos(|(vr + 2(^)) V ; 

where 5c^ k := Cj k — Y^°\ Of course, one has to keep in mind that these expressions 
are valid for small I. This corresponds to focusing on the kinematics near — ^13 — 
c i4 = c i4 i n the space of the cross-ratios (or equivalently near the regular decagon). 
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7. Comparison with two-loop results 

Wilson loops with light-like edges are dual to gluon scattering amplitudes [2]. For the 
kinematic configurations corresponding to the AdS$ octagon, the analytic expression 
of the 2-loop remainder function for the Wilson loop is given in [19]. The analytic 
expression for the case of AdS$ 2n-gon has also been written down [201121] • In this 
section, we compare our strong coupling results with those at two loops as expansions 
around the kinematic configurations associated with regular polygons. 



7.1. Octagon 

In the case of the octagon, the remainder function at two loops 



^2-ioop = - I log(l + x + ) log (l + log(l + X ~) log (l + ^r) 



(7-1) 



where the cross-ratios are given in (15. 6p . Similarly to the strong coupling case, 
this is expanded by using ( lA.lOj) as 



OC 



^2-loop = ^-1°°P (2*) (0) fk ; 



? 2-loop(2)^ _ log 2 (log 2 



fc=0 

where the first few coefficients are 

< loop(0) ( ^_log^2 

o 

pa-loop (4) 21og 3 2-51og 2 2 + 41og2-2 

^8 10J = I'-^J 

2 log 3 2 + 3 log 2 2 - 12 log2 + 6 

H cos 40. 

1536 v 

One can check that R% loop (°) agrees with i?g~ loop for the regular octagon [TH], and 
that the coefficients i?g loop ( 2fc ) maintain the space-time parity and cyclicity. 

For comparison of the results at strong coupling and at two loops, we introduce 
rescaled remainder functions [52]. For the AdS$ 2n-gon, they are defined by 

p . -R2n ~ R2n,reg (7 A \ 

2n -~ ~R _ ( n -9^R ' ^ ' 

at strong coupling, and a similar expression at two loops, where -R2n,reg stands for the 
remainder function for the regular 2n-gon. Since R 2n , R^n°° V reduce to superpositions 



The overall coupling dependence is suppressed. 
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of the contributions from (n — 2) regular hexagons in the low-temperature/collinear 
limit [HIED], the rescaled remainder functions are calibrated to take —1 in this limit. 
It has been observed numerically [H],[52] that R$ at strong coupling and -Rg loop at 
two loops are very similar. 

Given ( I5.15P and (j7.2j) . we are now able to derive analytic expansions of R$ and 
.Rg loop . By noting that the remainder functions for the regular hexagon and octagon 
are 

-R6,reg = TT) -R8,reg = R$ ^ ) 



and 

^)2-loop 

36 



4 

p2-loop _ _ 7F c>2-loop _ t-,2- 1o °p(0) (<-! a\ 

-"^reg — o/?' ""^reg ~~ -"-8 ' \'- v ) 



x: 



respectively, we obtain 

oo 

i? 8 = } (0)P, = ^ i? 2-loo P (2 fe)(0) p (? ?) 

fc=l fc=l 

where the first few coefficients at strong coupling are 
R^U) = — -if- « -0.1637687, 

8 VYV _ log 2 2 ' 

12 2 

i^ 4) (0) « 0.0174868 + 0.000667828 cos 40, (7.8) 
^8 6) 0) ~ -0.00160021 - 0.000173979 cos 40, 

whereas those at two loops are 

R% lo °vWU) = lo & 2 ~ 1 w -0.15966848, 
8 y ' 4 log 2 2 

J Rg" loop(4) (0) « 0.0163067 + 0.00118658 cos40, (7.9) 
J Rg" loop(6) (0) « -0.00141679 - 0.00029145 cos40. 

We observe that they are indeed close to each other (but different). 
7.2. Decagon 

Let us move on to a discussion on the AdS% decagon. In this case, the analytic 
expression of the remainder function in [2U1EI] is 

^2-loop = -— - - log(Ufc) log(ujH-l) log(M fe+2 ) log(w fc+3 ) , (7.10) 
k=l 



35 



with Uk = Wfc+io- The cross-ratios Uk are related to cf 3 , cf 4 b} 



«2 = 
«4 = 
M 6 = 
U8 = 



1 + C 



+ 
13 



1 + C+ +C+ ' 



1-1 



-11 



C++1 



1 + C+ + C 



11 



;i + 4)(i + c+) 



1.3 



C++1 



1 + C 



11 



l + cj 3 + c 



-11 



Ml 

^3 

U 7 
U 9 



1 + c 



13 



1 + c 13 + C 



11 



11 



C 14 + 1 



1 + c 13 + c 



11 



;i + cr 3 )(i + cr 4 ) 

c 13 



(7.H) 



c 13 + 1 



1 + c 



11 



1 + c 13 + c 



11 



Since i?i loop is invariant under the symmetries associated with the space-time 
parity and cyclicity, its high-temperature expansion is largely constrained similarly 
to AAbds- Consequently, by substituting the cross-ratios (16. 5p . (16. 6 p into (I7.10p . 
(17. lip , one obtains the following expansion of the remainder function at two loops: 



1-,2-loop _ V"^ p 2-loop (fc) ,2fc/5 



(7.12) 



fc=0 



where the first few coefficients are 



R 



R 
R 



2-loop (0) 
10 

2-loop (1) 
10 

2-loop (4) 
10 



-^-51og 4 (2cos^. 



R 



2-loop (2) 
10 



R 



2-loop (3) 
10 



(7.13) 



D 2 ■ \Y {2) {M x e iip \ M 2 e 



with 



2 4 x/5 



cos 



2n 



log I 2 cos 



7T 



3^ 



2 4 cos 2 



log I 2 cos 



7T 



5/J 



(7.14) 



One can check that i?^ loop ^ agrees with the numerical value of R^ oop for the regular 
decagon [52]. The structure of the expansion is very similar to that of Riq at strong 
coupling, which is understood as due to the space-time symmetries (see Appendix D 
for details.). 

Given the above result, we can compare the rescaled remainder functions at strong 
coupling and at two loops. From (I7.4p and a similar expression with 



R 



10,reg 



R 



(0) 
10 ) 



p 2-loop 
-"10,reg 



R 



2-loop (0) 
10 



(7.15) 



9 We identify in [20 with xjf , so that the Zio symmetry from the parity and cyclicity matches 
at strong coupling and at two loops. 
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Figure 4: Plots of the /-dependence of the rescaled remainder functions at strong 
coupling (points) and at two loops (dashed lines). The functions are evaluated at 
M 1 = M 2 = 1 and (p ± = cp 2 = tt/20. 



we find that 



R 10 = C 8/5 |y( 2) (Mie^M 2 e l ''<' 2 )| 2 ■ Z 8/5 + C(/ 12/5 ) , 
r^p = Cl'/° op \Y^\M 1 e itp \M 2 e i ^)\ 2 ■ l 8/5 + C(/ 12/5 ) , (7.16) 



where 



-itanf + fia 

C 8 = 1 s * w -0.0441916, 

« f-|log 2 (2cosf) 

C'«" l00P = c , 4~.f 2 Ti- \ ~ -0.0449039 . (7.17) 
5 5 log (2 cos ! J 

B2 and are given in f)6.30p and f)7.14p . respectively. Again, we observe that 
they are very close. We note that the two functions are also very close for finite I 
(Fig. H]). This suggests that not only the high-temperature expansion but also the 
remainder function itself is strongly constrained by the Y-system and the space-time 
symmetries, in addition to the collinear limits [20~ | l2~l~ | l53] . 

8. Conclusions and discussion 

In this paper we have studied the remainder functions of the gluon scattering ampli- 
tudes at strong coupling by using the integrable bulk and boundary perturbation of 
conformal field theory. In particular we have studied the minimal surfaces in AdS^, 
which correspond to the Wilson loops with a 2n-sided light-like polygonal boundary. 



37 



The minimal surface is described by the TBA system, and the related integrable 
model is the homogeneous sine-Gordon model with purely imaginary resonance pa- 
rameters. This model is obtained by the integrable perturbation of generalized 
parafermions. We have investigated high-temperature (small mass) expansion of the 
free energy, the T-functions and Y-functions of this model, which give the remainder 
function around the kinematic configurations associated with regular polygons. The 
high-temperature expansion of the free energy is calculated by the bulk perturbation 
of the CFT. For the T-functions, we have introduced the ^-functions whose ratios 
obey the same integral equations and asymptotic conditions. By using this relation, 
we have calculated the T-functions. 

For the 8-point amplitudes (octagon), the relevant CFT is the critical Ising model. 
Since we know the exact ^-function in this case, we have performed all-order high- 
temperature expansion of the remainder function. We have compared this result 
with the 2-loop remainder function, and have observed that the two results show 
similar power series expansions with very close coefficients. 

We have also been able to obtain an explicit formula for the first order correc- 
tion to the remainder function in the case of the 10-point amplitudes (decagon). 
The correction agrees with the numerical solution of the TBA equations for small 
masses. We have compared this result with the proposed 2-loop remainder function. 
Again, we have observed that the rescaled remainder functions have similar power 
series structures with close coefficients. This observed similarity suggests that their 
power series structure is strongly constrained by the Y-system and the space-time 
symmetries, in addition to the collinear limits [20 | [2T | l53] . 

The Y-functions are obtained from the cross-ratios of the T-functions. A notable 
observation in our discussion is that the T- and Y-functions in the CFT limit is 
given by the modular S-matrix through the ^-functions. This suggests an interesting 
relation between the modular S-matrices and solutions of constant Y-systems, which 
are used to compute the central charge of CFT using the dilogarithm identities 
[TUJSSlSn] - Although we have also observed that the same integral equations are 
derived from the g- and T-functions, the role of the boundary perturbation of CFT 
in the context of gluon scattering amplitudes is not quite clear at this moment. These 
points would deserve further investigations. 

For future direction, it would be possible to extend our discussion to the cases of 
more than 10-point amplitudes, amplitudes with more general kinematic configura- 
tions corresponding to AdS± and AdS&, and form factors [22J. For these purposes, it 
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would be important to understand multi-parameter integrable deformations of the 
generalized parafermionic CFT. In addition, the underlying integrable models/CFTs 
for the AdS*> case are yet to be clarified. Taking into account the case of the AdS§ 
hexagon [HII3], one may expect them to be obtained by some deformation from the 
AdSi case [9]. Details should, however, be discussed further. 

Regarding higher order expansions, the approach adopted in this paper requires 
higher correlation functions in the presence of both the bulk and the boundary de- 
formations. In order to cover the full kinematic region of gluon momenta, one needs 
higher order expansion, which is connected to the low-temperature (large mass) re- 
gion [7J|9j|53j|5l]. On the other hand, a different approach to study the analytic 
expansion of the T-functions has been given by Bazhanov et al. for kink (massless) 
TBA systems (see Appendix C) [3^|I53] . For massive systems, see for example [5B|I5T]. 
It would be interesting to apply this to the study of the minimal surface, as well as 
to understand the relation between these two approaches. 
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A. Expansion of F(l, cp) 

In this appendix, we derive the high-temperature expansion of F(l, <p) defined by 
( 15. 9p . First we expand F(l,<p) in sin (p, 
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COS if 



oo 



cosht 



F(l,<p) 




log(l + e 



—l cosh t 



7T 



cosh 2 1 — sin 2 ip 



oo 




(A.l) 



n=0 
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where 

MO 



dt 



-kl cosh t 



cosh m t 



(A.2) 



It is easy to see that the m-th derivative of $ m (i) is expressed in terms of the modified 
Bessel function of the second kind, 



$M(Z) = ^(-l) fc + m - 1 A; m - 1 ir (A;/). 



(A.3) 



k=l 



The summation (1A.3j) for m = 2n + 1 can be evaluated by using the following 
formula 



^(-l) fc cos(A;a)iro(^) 
k=i 

=K 7£+iog i) + ^ 



1 



1 



j=i L 



y/ 2 + [(2j - 1)tt - a] 2 2ttj 

1 1 

V/ 2 + P - 1)tt + a] 2 ~ 



For example, 



-Z 2 + 2(2j - 1) 2 tt 2 



,i= 

CO 



^ [/ 2 + (2j-l) 2 7T 2 ] 5 /2' 

3/ 4 - 24(2j - 1) 2 tt 2 / 2 + 8(2 j - 1) 4 tt 4 



3vrg [P + ( 2j - 1)2^9/2 

By integrating both sides in (IA.5j) . ( )A.6j) etc., we find the general structure 

2n <i) (j) rn^i 00 / 2 ™ 

«T,' fA — ^2n+ll u J , 7 -l / -,xn [ 



One can explicitly check this for small rt's. Therefore 

2 ™ 

J! 

1 



i=0 



m — n 



1 



22mH 



I C(2m+1) 



1 



2m + 1 V 7r 



2m+l 



(A.4) 



(A.5) 
(A.6) 



(A.7) 



(A.8) 
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The remaining task is to determine $2n+i(0) (j = 0, 1, . . . , 2n). If we define f(x) 
log(l + e~ x ), from (1A.2D we find 



*SUo) = / w (o) 



(3)1 



1 aJ /W(0). 



cosh 2n - J+1 t " r(2n-j' + l) 
From the following series expansion, 



. , x ^ (2* - 1)B 
log cosh - = ^ — — — 



2 ^ J!-J 



we get 



/(0) = log 2, /W(0) = (2 ' (j > 1), 



(A.9) 



(A.10) 



(All) 



where the Bernoulli numbers are defined by 



x 



B„ 



1 2s n\ 



n=0 



(A.12) 



Substituting (1A.9j) and ( lA.llj) into ( 1A.8|) . we obtain the expansion of $2n+i(0- I 11 
summary, the expansions of 3>2n+i(0 are 



7T 



$l(Z) = - log 2 - - ( logy + 1 + l0g7T - 7e 



2 

rn 



$2n+l(0 



1)2 
2^ 



m=l 
2 2n-l r ^ n + 

r(2n + 1) 

00 

+ (-1)"* £ 



2 2r 



I C(2m + 1) 



1 // 



2m + 1 V 7T 



2m+l 



(A.13) 



2h 



log2 + ^ 



2 2n-l( 1 _ 2-3) T( 



n 



J-1^2 



?B 3 



r(2n-j + l) 



V 



m — n 



I C(2m + 1) 



2m + 1 V 7r 



2m+l 



(A.14) 



Substituting flA.13[) and ( 1A.14|) into ( 1A.1[) . we obtain 



F(^) = ilog2-i- 



cos ( log - + 1 + log 7r - 7 E ) + (p sin 



p=2 

(A.15) 
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The expressions of Cp(ip) are given by 

(2 2k — l)B 2k . 2k 
CokW) = n rr. — sin (p, 

C 2k+1 ((f) 



= ^Vsin 2 >.(-irvrf ~i ) (l- 1 W + l)-^ i- 

7T ^ v ; \k-nj \ 2 2k + l ) sv y 2A; + 1 7r 2 * +1 

= (?) (^TTKTO (l - ^Tl) ^ + 1) - <p 2 F l( -k, 1; 1 - Jfc; sin 2 „). 

(A.16) 

Note that the sum over even p can be performed, 

£ c 2k ( V )l 2k = - ( I sin <p + 2 log -— . (A.17) 

fc=i ^ ' 

Thus we finally arrive at the expansion (15. lip . 

B. Perturbation of the (SU(2)_ 1/2 x SU(2) 1 )/SU(2) 1/a coset model 

In this appendix, we compute the high-temperature expansion of the free energy 
in the non-unitary (SU(2)_i/ 2 x SU(2) 1 )/SU(2) 1 / 2 coset model perturbed by the 
primary field 0i,i,3 with dimension A = A = 1/5. This model plays an important 
role in analyzing the remainder function for the decagon with Mi = M 2 . This model 
is equivalent to the non-unitary minimal model M.^,b perturbed by the relevant 
operator $ = $1,3, whose action takes the following form, 



S = S CFT + Xj d 2 x$(x). (B.l) 
Using the result in [SS], we can write down the coupling- mass relation 

A = kM 8/5 , 



k2 = - 2l [-5)\5 



7T ( 1 

7 <4 



1 16/5 

(B.2) 



Note that k is purely imaginary in this case as well as in the scaling Lee- Yang 
model Ai 2 ,5- The central charge of the UV CFT is c = — 3/5, and the ground state 
corresponds to the operator $i j2 with dimension Ao = Ao = —1/20. Thus the 
effective central charge is given by 

c = c-12(A + A ) = |, (B.3) 

5 
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which is one-half of the central charge for the SU(3) 2 /U(1) 2 coset CFT as expected. 

Let us consider the free energy of this model. Near the high-temperature limit 
I — y 0, the free energy is expanded as 



1 oo 

F{l) = lc+-/ + Y.f {n)l * n '*- (B.4) 

where 



6 4 

n=l 



f(n) = lfpc<- n \ (B.5) 
6 

-Tl n r r ^l rl 2 7- 



£(«) = 12 LiZV 27r )2A-i f TT 
n\ J J -_ J - 



1 = 1 (2^1^-1)2(1-^) 
x ($ (oo, oo)$(l, l)$(«i, zi) • • • $(z n _i, z n _i)$ (0, 0)) 

connected • 

(B.6) 

Recall that the vacuum operator is <3>o — $1,2 and the perturbing operator is $ = $1,3. 
The first non- vanishing coefficient is C^: 

= -12(27r)- 3 / 5 C Wo , (B.7) 

where C$ $$ is the structure constant, which was computed in [59J, 

(C m f = Pji'ala,*)) 2 = 7 (-\) 7 (§) 7 (|) 2 • (B.8) 
Thus the leading correction is 



' l "=i5U -UJ H-s; 7 UJ 7 UJ- (B ' 9) 

Returning to the discussion on the decagon, we obtain an expansion of the free energy 
with Mi = M 2 = M, 

At:r n \ Ml =M 2 = 2F(i) = \ + \l 2 + 2/ (1) ^ 8/5 + " " •■ (B.10) 
C. Generalization to complex masses 

Here we discuss how to extend the expansions of the Y-functions to the case with 
complex masses. We focus on the decagonal case n = 5. However the generalization 
to n > 6 is straightforward. 

Let us consider the expansion of the Y-functions. From the quasi-periodicity 
(14.211) and analyticity, the Y-functions should have the following expansion [11 J 



2 

k=— oo 
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with Y 2 = (—l) k Yj k \ In addition, the reality condition (12. lip constrains the 



coefficients as Yj ^ = Y-, which gives 

Y 1 (6) = \f2( Yik)ef0 + ¥me ~ f0 )> 

1 oo 

Y 2 (0) = i^(-l) fe (yW e fe + 7(fc)e-f e ) . (C.2) 



fc=0 

The coefficients Y^ are functions oirrij and fhj, and Y^ here coincides with (I6.25I) 
if all the masses are real. Y^ are expanded in powers of I 2 / 5 , with the leading 
behavior Y^ ~ / 2fc / 5 for small I [TTIHB]. Moreover, according to (10], this leading 
behavior is thought of as coming from the boundary perturbation in ( I4.16p . The 
form of the expansion ( I4.16P then implies that the subleading corrections are given 
in powers of / 4 / 5 . Thus, one may have 

oo 

>^ = $>, 2p /^+ 2p \ (C.3) 

p=0 

where I = ML and M is an overall mass scale. This is in accord with the double 
expansion in terms of le ±e discussed in jlH] . At low orders, the absence of terms of 
order / 2 ( fe +p')/5 w ith odd p' is also confirmed from the Y-system by following [31]. We 
have checked that the above expansion is consistent with numerical results. 

It is convenient here to write the first few terms of the expansion of Yi{6) in l 2 ^ 5 , 

2Y 1 (9) = 2b 00 + (b 20 ef + b 20 e-f)l^ + 0(l e / 5 ), (C.4) 

where we have used the fact that b 00 = b 00 = 2cos(7r/5) and b w = b 02 = 6 i2 = 
as seen in Appendix D. The coefficients b^^p depend on both Mje~ llfj and Mje liPj in 
general. It is important to notice that from (I2.17P the TBA equations for Yj(6) = 
Yj(8 — iipj) are given by 

log Y,(0) = - l -{m 3 e e + m 3 e- e ) + K * log(l + + Y j+1 ). (C.5) 

In order to reveal the complex mass dependence of bk,2n, we consider the decoupling 
limit (chiral limit) I -> 0. In this limit, the new functions Y J kmk (6 ) ) = Yj(8 — log(Z/2)) 
satisfy the kink TBA equations (see [12] for example) 

log Y^\6) = -Mje-^e 9 + K * log(l + Y*f )(1 + Y^f). (C.6) 
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From the periodicity, Y kmk (6) have the expansion, 

1 oo 

Y^\9) = -Y,y J kink{k)e2Jf - ( c - 7 ) 



2 

k=0 



(IC.6j) suggests that the coefficients y^ kmk ( fc ) are functions of Mje lVj , not of Mje llfj . 
On the other hand, by taking the decoupling limit in ( 1C.2j) with (1C.3j) . we obtain 

Y^\6) = limn (fl - log ^) = 600 + \ V 2" b k0 e 2J f. (C.8) 
™ V V 2^ 



Comparing (IC8[) with (ICT[) . we obtain 



ykink(O) = 26q0) y kink(fc) = 2 2fc/5 6fco _ (c g) 

These relations show that the coefficients bko (k > 1) depend on Mje~ %v \ not on 
MjC lip3 . Similarly, bko (k > 1) are functions of Mje 1 ^ . In summary, we can get 
bw (pko) for the complex masses by replacing Mj — > Mje~ t!P: > (Mje l{pi ) in bko for the 
real masses assuming the analyticity in Mj. However, b^^n { n > 1) are functions of 
Mje~ lLP: > and Mje l<Pj , and the above argument does not applied to them. We already 
know the small I expansion of the Y-functions for the real masses up to order / 4 / 5 
(see (ESI). The coefficient of is 

6 rcal = Y ^\M U M 2 ). (CIO) 

Using this result, we can obtain 6 2 o and 620 for the complex masses by the above 
prescription, 

^complex = Y ( 2 )(Mie- i(fi \ M 2 e^ 2 ), (C.ll) 

^complex = y^^Mxe^^M^ 2 ). (C.12) 

Substituting these equations into (IC.4I) . we obtain the small I expansion of the Y- 
functions for the complex masses as in (16.401) . 

D. Structure of expansions at higher orders 

In the main text, we have obtained the first order high-temperature expansion of the 
Y-functions for the decagon by using the conformal perturbation. In this appendix, 
we show that the structure of the high-temperature expansion is largely constrained 
by the Y-system, the structure of the conformal perturbation, and the symmetries 
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associated with the space-time parity and cyclicity, although one still needs higher 
order perturbations to find precise values of the coefficients. In the following, we con- 
centrate on the case of the AdS^ decagon, but the discussion below can be extended 
to more general cases. 

We start with the expansion (IC.2j) . and flC.31) which is in accord with the confor- 
mal perturbation as discussed in Appendix C. Substituting these into the Y-system 
(12.101) for the AdS3 decagon, one obtains a double expansion in e 2e ^ 5 and Z 2//5 , in 
which each coefficient should vanish. For the first few orders, we then find, e.g., 

600 = 2 cos(-) , 610 = 6 02 = 612 = , 

6 04 = - sin 2 (-^) ■ 620^20, ho = jjsin 2 (^) • (6 20 ) 2 • (D.l) 

The expansion of the Y-functions in turn gives the expansions of the cross-ratios and 
AAbds- Using the relations among bk )2p obtained from the Y-system, we find that 



AAbds = Yl A ^ ' (D.2) 



k=0 



with 



A = -^log 2 (2cos^), A 4 = B 2 ■ b 20 b 20 , (D.3) 

2 

A 6 = 5 2 (620622 + 622620 + 630630) - 40 cos 4 (-71-) ■ 6 30 6 3 o , 

and A\ = A 2 = A 3 = A 5 = A 7 = 0. B 2 is given in (16.301) . The result shows that one 
can obtain the expansion of AAbds u P to (9(/ 12//5 ) once 6 2 o is known. In the main 
text and Appendix C, 6 20 is found to be 6 20 = F (2) (Mie - ^ 1 , M 2 e"^ 2 ) . 

One can argue that the absence of A& with odd k is understood as a consequence 
of the Z10 symmetry due to the space-time parity and cyclicity, x^ — > xf +1 , xf — > 
The Z10 transformation is concisely expressed by the Y-functions as [21] 



Yj (9)^ Yj (9+^i), (D.4) 

or in terms of the expansion coefficients and the TBA masses, — > e^ k Y^ 
and rrij — > rrij/i, respectively. The cross-ratios cf 3 ,cf A transform as c^ 3 — > c^, 
c u ~^ c 25' c i3 ~^ C ni c i4 c ui where c^jC^ are given in (16 77$ . Av4 B ds in 
(16. 4p indeed has this symmetry. Note also that A^ consists of terms of the form 
II V,,2 gj with k = J2(pj + 2qj) + J2(Pj + 2 ?i); which transform under f lD.41) 

as II b Pj ,2 qj II bpjMj -> e ^ Ap 11 b Pj ,2 qj ll h,M, with A P '■ = E Pi - E Pj- Thus, unless 
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non-trivial cancellations with other terms occur, which is unlikely because are 
functions of rrij,rhj, the terms with Ap ^ (mod 10) are projected out by successive 
actions of the Zio transformation. For odd k, Ap necessary becomes non-zero and 
hence is projected out. This leads to the symmetry / 2//5 — > — Z 2 / 5 . Furthermore, 
the Zio symmetry is then promoted to a continuous symmetry — > e t( t> k Y^ k \ 
which in turn corresponds to a continuous imaginary shift of 9 in Yj(9) and hence 
the world-sheet rotational symmetry. 

AAbds is also invariant under the space-time parity symmetry: xf — > x~ with 
the order of the labeling of the cusps being reversed. In the case of the decagon, 
this gives rise to c^ 3 -H- c u , which is equivalent to l 2 ^ — > — Z 2 / 5 and Y^ -H- Y^ k \ 
Combined with the above symmetry Z 2 / 5 — > — Z 2//5 , the parity results in the symmetry 
^ Consequently, only the terms with Ap = and maintaining the 

symmetry as b P: 2qb P 2q> + b P: 2qb P: 2q', are allowed in the expansion. 

In addition, since the remainder function has to have the parity and Z i0 symme- 
tries, v4 period + A {ree (= R w — Av4 B ds) also maintains these symmetries. Oppositely, 
since A period + A ivee is invariant under the world-sheet rotational symmetry by defi- 
nition, so is the remainder function. 

Regarding the functional form of b p ^ q , the conformal perturbation and the expan- 
sion in le ±e [46J for complex TBA masses suggest that b p ^ q are given by summation 
of terms of the from l~^ p+2q ^rh 2 / 5 ■ ■ ■ fh 2 / 5 ■ m 2 ! 5 ■ • ■ m 2 J 5 . Indeed, given this form, one 

31 3p+q 3x 3q ' ° 

finds that Y^ transform as Y^ — > e^ k Y^ and Y^ — > e~ % ^" p Y ( ^ 1 under the Z 10 
transformation rrij — > rrij/i and the world-sheet rotation equivalent to rrij — > e llp mj, 
respectively. This is in accord with (ID.4[) and the above argument. Note that the 
invariants under the latter symmetry depend on the phases of rrij only through their 
deference </?i2 = ipi — if2- 

Our results of the expansion are all consistent with the above arguments. How- 
ever, further investigations are needed for definite conclusions. 

We can also study the structure of the high-temperature expansion of the remain- 
der function at two loops. From the expansion of Yj(9) in (1C.2|) . (1C.3I) . we find that 
the remainder function (I7.10P is expanded as 



jR 2- ) loop = ^ i? 2-loop( fc )p / 5 ; (D5) 
fc=0 
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with 

iC° p(0) = ~ - 5 log 4 (2 cos |) , iC ^ = D 2 ■ b 20 b 20 , (D.6) 

^2-loop (6) = D ^ bw l 22 + ^ + hQ - hQ) + 5 . 2 8 CQS 9 ^ lQg 2 ^ CQS ^ . ^ ; 

and i4 Ioop(1) = it!j loop(2) = itf loop(3) = itf loop(5) = itf loop(7) = 0. D 2 is given in 
(I7.14p . Thus, the structure of the expansion is very similar to that of AAbds- This 
is because i?^ loop has the parity and Zio symmetries and hence the structure of the 
expansion is strongly constrained as in the case of AAbds- 
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